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Abstract 


pEst parameter estimation (NASA Dryden code) 


The application of an adaptive real-time measurement- PSC 
based performance optimization technique is being 
explored for a future flight research program. The key Symbols 
technical challenge of the approach is parameter identifica- 
tion, which uses a perturbation-search technique to iden- a 
tify changes in performance caused by forced oscillations 
of the controls. The controls on the NASA F- 15 highly AI T 

integrated digital electronic control (HIDEQ aircraft were 
perturbed using inlet cowl rotation steps at various sub- u 
sonic and supersonic flight conditions to determine the 
effect on aircraft performance. The feasibility of the 
perturbation-search technique for identifying integrated c 
airframe-propulsion system performance effects was sue- L 
cessfully shown through flight experiments and postflight C 
data analysis. Aircraft response and control data were ana- " 
lyzed postflight to identify gradients and to determine the 
minimum drag point. Changes in longitudinal acceleration ° 
as small as 0.004 g were measured, and absolute resolution C m„ 
was estimated to be 0.002 g or approximately 50 lbf of 
drag. Two techniques for identifying performance gradi- 
ents were compared: a least-squares estimation algorithm ^ m s 
and a modified maximum likelihood estimator algorithm. 

A complementary filter algorithm was used with the least 
squares estimator. 


performance seeking control 


estimated coefficient of output series in least- 
squares analysis 

pressure altitude, ft 


estimated coefficient of input series in least- 
squares analysis 

total lift coefficient 

total pitching moment coefficient 

pitching moment coefficient bias 

derivative of pitching moment coefficient 
with respect to angle of attack, deg -1 

derivative of pitching moment coefficient 
with respect to symmetric cowl deflec- 
tion, deg -1 


Nomenclature 


derivative of pitching moment coefficient 
with respect to stabilator deflection, 


Acronyms 

DEEC 

digital electronic engine control 

S 

— O 

derivative of pitching moment coefficient 
with respect to pitch rate, rad -1 

HIDEC 

highly integrated digital electronic control 


total longitudinal force coefficient 

LSE 

least-squares estimator 

C *0 

longitudinal force coefficient bias 

MLE 

maximum likelihood estimator 

c *. 

derivative of longitudinal force coefficient 
with respect to angle of attack, deg -1 
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derivative of longitudinal force coefficient 
with respect to symmetric cowl deflec- 
tion, deg -1 
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derivative of longitudinal force coefficient 
with respect to stabilator deflection, 

deg -1 

total normal force coefficient 


dynamic pressure, lbf/ft 2 
wing area, 608 ft 2 
Laplace operator 


normal force coefficient bias 


derivative of normal force coefficient with 
respect to angle of attack, deg 

derivative of normal force coefficient with 
respect to symmetric cowl deflection, 
deg -1 

derivative of normal force coefficient with 
respect to stabilator deflection, deg 

noise term in least-squares analysis 


excess thrust, lbf 


MLE state function 
transfer function from u to y 


high-pass component of complementary fil- 


ter 


low-pass component of complementary filter 

gravitational acceleration, 32.174 ft/sec , 

6 also function for MLE output 


transfer function from e to y 

likelihood cost function 

Mach number 

mass of the aircraft, slugs, or pitching 
moment 

longitudinal acceleration, g 

high-pass filtered body-axis longitudinal 
acceleration, ft/sec 2 

complementary filtered longitudinal acceler- 
ation, ft/sec 2 

normal acceleration, g 

number of delay terms in input-output rela- 
tionship 

pitch rate, deg/sec 

delay operator 
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time 

input parameter in LSE analysis, also control 
vector in MLE analysis 

velocity, ft/sec 

derivative of low-pass filtered body-axis 
velocity, ft/sec 2 

derivative of longitudinal force with respect 
to symmetric cowl deflection, gl deg 

state vector in MLE analysis 

output parameter in LSE analysis, also com- 
puted observation vector in MLE analysis 

one-step-ahead prediction of output parame- 
ter in LSE analysis 

measured observation vector in MLE 
analysis 

angle of attack, deg 

perturbation of control 

control surface deflection, deg 

damping ratio, 0.707, also parameter vector 

noise vector in MLE analysis 

pitch attitude, deg, also vector of estimated 
coefficients in LSE analysis 


<p transpose of input and output data vector 

w natural frequency, 0.1 rad/sec 

Subscripts 

a last output coefficient in LSE input-output 

relationship 

b last input coefficient in LSE input-output 

relationship 


complementary filtered 

symmetric cowl deflection, positive leading 
edge down, deg 
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e symmetric stabilator deflection, positive 

trailing edge down, deg 

i index of discrete cost function in MLE 

inf free-stream static 

k number of delays from input to output in 

LSE input-output relationship 

m measured 

0 trim position 

Superscripts 

T transpose 

Introduction 

For over a decade, the NASA Dryden Flight Research 
Facility has conducted a multidisciplinary performance 
improvement flight research program on an F- 15 aircraft. 
Significant parts of this research involved the flight evalu- 
ation of advanced propulsion control concepts in pro- 
grams such as digital electronic engine control (DEEC), 1,2 
the F100 engine model derivative (EMD), 3 and highly 
integrated digital electronic control (HIDEC). 4 - 5 The 
increased performance and improved fuel economy shown 
on the F-15 HIDEC research vehicle led to the perfor- 
mance seeking control (PSC) program in which adaptive 
optimization of the near steady-state performance of an 
aircraft-propulsion system was addressed. 

The PSC approach computed and applied incremental 
trims to the aircraft-propulsion control system input com- 
mand schedule to enhance the performance of the air- 
craft. 6 ’ 7 The PSC control law included an estimation 
process using a Kalman filter to update the onboard engine 
model to represent the actual engine hardware better. The 
optimization process used linear-programming techniques 
to determine the optimal engine operating condition for 
the selected performance measure (Fig. 1). Although the 
PSC algorithm showed substantial benefits, there were 
inherent shortcomings. The PSC approach to optimizing 
aircraft performance consisted of a model-based, open- 
loop algorithm. The algorithm relied heavily on accurate 
models of airframe, engine, inlet, and nozzle for the opti- 
mization process. No intrinsic means in the PSC approach 
could compensate for errors in modeling and measure- 
ments. The Kalman-filter estimator depended on accurate 
engine measurements and could not account for biases in 
those measurements. 8 

Extensive testing of the PSC algorithm for the subsonic 
phase and preliminary testing of the supersonic phase 
have already been successfully conducted. 9 * 10 ’ 11 Right 
test results have validated predicted improvements in fuel 


efficiency, increases in thrust, and reductions in fan turbine 
inlet temperature leading to longer engine life. 

The subsonic phase of the PSC program dealt exclu- 
sively with engine optimization, while the supersonic 
phase included inlet and stabilator integration. Optimiza- 
tion was based on models of the engine, nozzle, and inlet 
and a predicted minimum drag relationship between the 
inlet cowl and the stabilator. The adaptive features of PSC 
dealt only with the engine model and did not address 
updating of inlet or airframe model characteristics. An 
altitude control outer loop indirectly controlled the 
stabilator. 

To address many technical shortcomings of the PSC 
approach, NASA Dryden is studying augmenting the real- 
time optimization problem with flight measurements and 
feedback control. 12 The necessity for costly, high-fidelity 
engine and aircraft models required for the PSC program 
has been lessened by using a direct measurement approach 
that determines the sensitivity of the performance 
parameter — i.e. t excess thrust or range factor — to control 
perturbations (aircraft and propulsion). Figure 2 depicts 
the direct measurement approach applied to drag minim i- 
zation on the F-15 aircraft Total aircraft drag is plotted as 
a function of stabilator position and cowl position for a 
given cruise condition. The scheduled operating point of 
the cowl dictates a certain stabilator position to maintain 
trim at the cruise condition. If a perturbation-search tech- 
nique is applied to the cowl to identify the performance 
gradient, the cowl can be trimmed, forcing the stabilator to 
a new trim position, resulting in a lower total drag configu- 
ration (higher excess thrust). 

The advanced airframe-propulsion technology control 
system concept, which uses this direct measurement 
approach, is divided into four major components: an exec- 
utive module, excitation module, parameter identification 
module, and controller (Fig. 3). The executive provides 
most importantly the fault detection logic such as channel 
comparisons and output limit checking. The executive 
also serves as the distribution point for all needed data to 
the other three modules. The excitation module sizes the 
perturbation signals — i.e., duration, amplitude, and 
frequency — based on information received from the exec- 
utive. The parameter identification module forms the par- 
tial derivatives of the performance parameter with respect 
to the controls. The controller receives the partial deriva- 
tives of the performance parameter with respect to the con- 
trols and applies trims to the controls in the favorable 
direction of reduced drag and increased fuel economy. 

Two major technical issues with this concept are its abil- 
ity to identify the required performance sensitivity with 
explicit control excitation that is not objectionable or 
noticeable to the pilot and whether conventional sensors 
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can measure the resulting effects. These abilities are con- 
sidered necessary for a real-time adaptive system. As 
such, a postflight evaluation of these issues was performed 
at selected regions of the F- 15 subsonic and supersonic 
envelope. 

This paper presents the feasibility of adaptively measur- 
ing in-flight performance sensitivities based on identifying 
changes in excess thrust with inlet cowl position. Data for 
this study were generated from flight test of the NASA 
F-15 HIDEC aircraft. Two algorithms are compared for 
identifying performance gradients: a least-squares estima- 
tion (LSE) and a maximum likelihood estimation (MLE) 
algorithm. The paper also describes the use of a comple- 
mentary filter algorithm to enhance the acceleration mea- 
surement for the LSE. 

Aircraft Description 

Performance sensitivity was studied on the NASA 
F-15A HIDEC research aircraft, S/N 835 (Fig. 4), which 
is a high-performance military fighter aircraft capable of 
speeds faster than Mach 2. Two F100 derivative 
(PW1128) afterburning turbofan engines power the NASA 
F-15 aircraft. The aircraft was modified with a digital 
electronic flight control system. More information on the 
F-15 can be found in Ref. 1. 

The PW1128 engine is a low-bypass ratio, twin-spool, 
afterburning turbofan technology demonstrator, derived 
from the F100-PW-100 engine. A full-authority DEEC 
similar to the one for current production F100-PW-220 
engine controls the engines. The DEEC software has been 
modified to accommodate PSC trim commands, but the 
normal DEEC control loops (i.e., corrected fan speed and 
engine pressure ratio) were not modified. The DEEC trim 
commands for subsonic, nonafterbuming conditions are 
perturbations on compressor inlet variable vanes, rear 
compressor variable vanes, corrected fan speed, and noz- 
zle throat area. A more detailed description of the 
PW1128 engine can be found in Ref. 3. 

The NASA F-15 aircraft was also modified with an 
electronic air inlet controller that allows PSC trim com- 
mands to be added with first and third inlet ramp sched- 
uled positions (Fig. 5). These inlet ramp schedules were 
tailored specifically for the PW1 128 engines during super- 
sonic flight to account for the higher engine airflow. 

Longitudinal acceleration (N x ) data were gathered from 
two longitudinal sensors, both flight test instrumentation. 
One accelerometer was in the noseboom, the other was 
near the center of gravity of the aircraft. The engineering 
unit ranges and resolutions of the accelerometers were 
±1.37 and ±1.03 g and 0.00268 and 0.00201 gl bit, 
respectively, using the aircraft's 10-bit digital-to-analog 


instrumentation system. The aircraft also had the standard 
set of stability and control and airdata parameters mea- 
sured. The data were recorded at 40 samples/sec, filtered 
and sampled down to 10 and 4 samples/sec for the post- 
flight analysis. 

Analysis Method 

Two approaches were evaluated to identify the sensitiv- 
ity of excess thrust with respect to changes in cowl posi- 
tion. The first method was a more direct approach, 
incorporating a one-degree -of- freedom , single-input- 
single-output LSE. To approximate the three-degree-of- 
freedom problem as having only one degree of freedom, 
the following external constraints on pitching moment and 
lift coefficients were imposed: 


and 

C L = constant 

The implementation of these constraints is described in a 
later section. Flight Test Maneuver Description and Condi- 
tions. 

The investigation assumed that changing the cowl posi- 
tion generates measurable performance changes. These 
performance changes result from a combination of aerody- 
namic forces and changes in engine thrust This coupling 
defines the integrated nature of the problem. The inlet 
geometry dictates the pressure recovery seen by the 
engine, which strongly affects thrust This coupling 
makes the problem difficult to model. 

In trim, the total pitching moment coefficient C m , on 
the aircraft is zero. 

C» = C «,' fC M a a + C m,« + C ”, > 8 . 

+ C -. 5 « = ° 

CM 

where 

C = pitching moment coefficient bias 

C = derivative of pitching moment coefficient 
m “ with respect to angle of attack 

C — derivative of pitching moment coefficient 
with respect to pitch rate 

C = derivative of pitching moment coefficient 
m *« with respect to stabilator deflection 
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^m 6 - derivative of pitching moment coefficient 

" with respect to symmetric cowl deflection 

If the aircraft configuration and angle of attack are 
assumed to remain constant, then changes (A) in cowl 
deflection will cause a stabilator trim position change to 
maintain the total pitching moment coefficient equal to 
zero. 


A 8, = -AS 
* CSC 


( 2 ) 


This new steady-state stabilator and cowl position is 
reflected in the longitudinal force coefficient, C x , eq uatio n 


“(0 = 8 cj , symmetric cowl deflection (first ramp) 

The Appendix develops the autoregressive exogenous 
model used for the LSE analysis. This model was chosen 
because it is a fast and uncomplicated algorithm that does 
not require an iterative solution. These characteristics 
allow the algorithm to be easily applied to a real-time sys- 
tem. Equation (5) shows the drag model used in the LSE 
analysis. The model consists of one output, three inputs (a 
constant, linear, and second-order term in u (/)), and no 
delays. 

Once the coefficients are estimated, the linear dimen- 
sional derivative, with units of gravity per degree, can be 
obtained about the trim condition by differentiating equa- 
tion (5) with respect to u (r) 


where 



( 3 ) 


dX 


06 


dyU) 

du(t) 


= b \ + lb 2 u Q 


( 6 ) 


where Uq = trim cowl position. 


C x 0 = longitudinal force coefficient bias 

C, Q = derivative of longitudinal force coefficient 
with respect to angle of attack 

c z h ~ derivative of longitudinal force coefficient 
* with respect to stabilator deflection 

C Xi = derivative of longitudinal force coefficient 
" with respect to symmetric cowl deflection 

and is directly measured by the longitudinal accelerometer 


N. 



Least Squares Estimator 


(4) 


The linear, time-invariant system model was made sec- 
ond order with respect to the first ramp position because of 
the second-order nature of drag with respect to wing inci- 
dence angle. No dynamic elements were introduced into 
this model, since forces imparted to and by the cowl were 
assumed to result instantaneously in acceleration. The 
output parameter, y (t ) , was chosen to be longitudinal 
acceleration. 

y(0 = b Q + b x u(t) +b 2 u 2 (t ) +e(t) (5) 


Complementary Filter 

To enhance the observability of the output parameter in 
the LSE analysis, a complementary-filter algorithm was 
used during the postflight processing of this investiga- 
tion. Complementary filters blend different sensor outputs 
to more accurately estimate the desired signal. For this 
investigation, corrected free- stream velocity, V- n ^, an air- 
data parameter, and longitudinal acceleration, N x , sensed 
by a force-balance accelerometer, were combined to pro- 
duce a complementary -filtered longitudinal acceleration, 
N *cf • with *e high-frequency characteristics of the accel- 
eration sensor and the low-frequency characteristics of the 
rate sensor (Fig. 6). This complementary filter in effect 
eliminated the bias normally seen with conventional accel- 
erometer data but retained the useful, high-frequency char- 
acteristics of the accelerometer. 

The complementary filter used for this study was made 
second order. The philosophy behind the complementary 
filter is that numerator and denominator polynomial 
coefficients are identical, so the magnitude and phase char- 
acteristics remain unchanged. 

n , . J 2 + 2£(oj + g) 2 

G(S) = “TT? j = 1 ( 7 > 

S + 2^0)5 + GT 


where 

y (0 — complementary longitudinal accelera- 
tion 


This polynomial is rewritten as the 
als. 


sum of two polynomi- 
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G(s) = G, (s) + G 2 (j) 


( 8 ) 



where 


G, (J) = 


S (5 + 2fo) 
s 2 + 2^(05 + (0 2 


(9) 


and 


0) 


G 2 (5) + + (0 2 


( 10 ) 


By examining these polynomials, one sees that G x {s) 
resembles a high-pass filter and G 2 (s) ™«***^ s of 
pass filter. The sensor characteristics of the two so 

used for the complementary filter match weUwth 

attributes of the filter. The force balance accelerometer 

has a linear, high-frequency response but « suba«t te 
hiases and low-frequency temperature drifts. The airdata 
SSjS toed^on measured static and total pressures 
resolve veJy small changes but is subject to lag 
and attenuation when high-frequency pressure c ang 
occur. It is natural, therefore, to pair the longitudinal 
accelerometer with G x (s) and the body-axis vetocity w 
Q (,) To obtain compatible units with acceleration, the 
derivative of the body-axis velocity must be taken; this is 
easily done by multiplying G 2 (s) by s, so G 2 (s) 
becomes 


FiKure 7 contains a block diagram of the program. Contra 
inputs are sent to the aircraft model generating an esu- 
mated response. The estimated response of the ancraft 
model Subtracted from the measured response * pro- 
duce the response error, which is sent to the MLE. The 
MLE has three components: a Gauss-Newton computa- 
tional algorithm, cost function, and parameter estimates. 

The Gauss-Newton computational algonthm is used to 
find the coefficient values that maximize the likelihood 
cost function. Each iteration of this algorithm provides 
revised estimates of the unknown coefficients 
response error. The revised estimates of the coefficients 
aredien used to update the aircraft model, providing a 
revised estimate response and, therefore a revised 
resoonse error. The mathematical model is updated ltera 
hvX until a convergence criterion is satisfied. For a more 
detailed discussion of the MLE approach and pEst see 
Refs. 14 to 17. 

A simplified MLE approach contained a three-degrees- 
of-freedom model that can estimate all long* jdinal nor idi- 
mensional derivatives or coefficients. The longitudinal 
aircraft equations can be written as follows. 


SCO 


G, (s) — “ 2 

2 s 2 + 2^(DJ+te‘ i 


( 11 ) 


id) =f[x(t),u(t),%] 
y(t) = g [x (0 . n (0 * 
z(r,) = >(*,-) +T t( r i) 


( 12 ) 


The frequency and damping ratio of the complementary 
filter were selected to be 0.1 rad/sec and 0.707, respec- 
tively, based on a brief study of the MDEC instrumenta- 
tion system. 

Maximum Likelihood Estimator 

The second method used in the study for comparison 
with the LSE method was an MLE. The MLE is a 
mathematically more precise probabilistic formulation 
than the LSE. For each possible estimate of the unknown 
parameters, a probability that the computed airc t 
response time histories attain values near the measured 
values can be calculated. When the estimates are chosen to 
maximize the probability of occurring, this is the maxi- 
mum likelihood estimate. Under certain conditions the 
MLE has many desirable characteristics, such as yi g 
asymptotically unbiased and consistent estimates. If the 
measurement noise is assumed to be Gaussian, white, «a- 
tionary. and uncorrelated, this formulation is equivalen to 
die LSE, where the weightings used are the inverse of the 
measurement noise covariance matrix. The MLE program 
used for the analysis was the NASA Dryden code, param- 
eter estimation (pEst). 


where 

x(t) - state vector 
M (f) = control vector 
y (j) = response vector 

z (t.) = measured output vector 

r\ (r-) = noise vector 

£ = parameter vector 

In vector notation the state, control, and observation vec- 
tors can be written as follows: 
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(13) 


The cost function can be defined as the integral squared 
error criterion and approximated in the discrete case as 
follows: 


6 



/( 5 ) = 


1 

2n z n , 


*1 

X !"(*(',) -y(ij) T 

1-1 L 


W [*(',•> ->(/,)] 


(14) 


where n t and n z are the number of time history points and 
response variables, respectively; t is the time variable; W 
is the response weighting matrix; 2 is the measured 
response; y is the response computed by integrating the 
equations of motion; ^ is the parameter vector; and super- 
script T denotes transpose. To obtain maximum likelihood 
estimates for a given set of flight responses, the cost func- 
tion, /, must be minimized. Since 2 (/) is fixed for given 
flight responses, J must be minimized by selecting the y 
that minimizes / where y is subject to the constraints of 
eq- (12). 


The nondimensional aerodynamic model used for the 
analysis was the standard three-degrees-of-freedom model 
contained in the pEst program's equations of motion. 


C + C ct + C 8 + C 8 

*0 x a x 6 * x t CS 

(15) 

* ct 


C 2 + C. a + C 8 + C 8 

z 0 z a *6 * *6 cs 

(16) 

* ct 


C m a + CV. <7+C 5 

m 0 m a m i e 


+ C 8 

(17) 

m b cs 

ct 



Although it would have been easy to augment the model 
by introducing second-order terms in controls, linear 
terms were sufficient and valid about the trim condition. 

In practice, not all coefficients were estimated at once. At 
each flight condition, a conventional doublet was first ana- 
lyzed to establish the basic stability and control coeffi- 
cients or derivatives: 
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These derivatives were then fixed during the analysis of 
the cowl step maneuver, and the remaining derivatives— 
the cowl control powers and biases — were estimated. 

This two-step approach prevented the dominant deriva- 
tives from trading with the cowl sensitivities. 

Flight Test Maneuver Description and Conditions 

To extract longitudinal acceleration gradients with 
respect to cowl deflections in flight required a well- 


conditioned maneuver. External perturbations were mini- 
mized to achieve the most accurate identification. As 
mentioned in the Analysis Method section, two external 
constraints were maintained during the maneuver: 


c m = o 

and 


C L = constant 

For a given Mach number ( M) and constant altitude 
(ALT), both of the above equations were constrained. The 
constant altitude actually was obtained using the altitude- 
hold autopilot mode when possible. The pilot selected the 
proper flight control mode, inlet integration, and trim val- 
ues to be sent to both cowls simultaneously. Next, the 
pilot stabilized the aircraft on condition in a hands-off, 1-g 
wings-level trim and engaged the autopilot At this point, 
stabilized trim data were gathered; no other aircraft config- 
uration changes were allowed. After approximately 30 sec 
of stable data, the pilot would select the positive cowl trim 
for 30 sec. Afterward, the pilot deselected the positive 
trim and repeated the procedure with the negative trim. 

Once the cowl step was completed a conventional longi- 
tudinal doublet usually was performed at the same condi- 
tions. The longitudinal doublet was used to estimate the 
lift and pitching moment derivatives of the aircraft at the 
conditions where the cowl step was performed. For the 
MLE analysis, these derivatives were held fixed to ensure 
that no residual forces entered the estimate. Figure 8 sum- 
marizes the test conditions flown for the investigation. 

Only three points were flown; however, they span the 
region of interest 

Results and Discussion 

The cowl was perturbed on the NASA F-15 HIDEC air- 
craft at three flight conditions to determine sensitivities of 
excess thrust to cowl deflections. Flight data were ana- 
lyzed postflight using two estimation algorithms for com- 
parison purposes. 

Figure 9 presents the results of a maneuver that exhib- 
ited the greatest change in excess thrust from the nominal 
schedule at a flight condition of Mach 1.25 at 25,000 ft 
Time histones are presented for measured performance 
and state parameters (excess thrust (F ex ), M , ALT, a, 0, q, 
N x , and N z ) and control parameters ( 8 cj> 8 ) . Stabilized 
trim data were gathered for the first 20 sec. The steady- 
state value of measured N x was approximately 0.01 g indi- 
caung a bias in the sensor or that the aircraft was in a 
slightly accelerating condition before the perturbations. 

The nominal cowl deflection and stabilator trim position 
were 4.7° and 4.3°, respectively. From 20 to 48 sec, the 
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first trim of -3° was sent to the cowl, causing the autopilot 
to drive the stabilalor 0.6° trading edge down to a new 
trim position of 4.9°. This reduced N x and F tx by 0.008 g 
and 500 lbf, respectively, and increased the trim drag on 
the aircraft. From 48 to 70 sec, the trim was removed and 
a dditional stabilized data were gathered. From 70 to 90 
sec, the last trim of +3° was sent to the cowl, causing the 
autopilot to drive the stabilator trailing edge up 0.55° to a 
new trim position of 3.85°. This trim increased N x and 
F by 0.023 g and 800 lbf, respectively, and decreased 

the trim drag on the aircraft. 

Since ride quality will limit the size of the disturbances 
from the desired flightpath, it is important to quantify the 
effects on ALT. 0,0, <?, and N,. Excursions of 14 ft tn 
ALT, 0.17° in a, 0.48° in 9 . 0.3°/sec in q, and 0.05 g in 
N occurred because of the perturbation. The pilot did not 
detect the cowl excitation and its effects. Thus the pertur- 
bation required to implement an adaptive control system 
of the cowl would not distract the pilot or passengers. 

Figure 10 shows plotted time histories of measured N x ^ 
and those computed using the estimated coefficients of the 
LSE algorithm. The computed and measured time 
histories match well with no bias or initial slope errors. In 
Fig. 11, time histories are plotted from measured flight 
data and computations made by the MLE algorithm (same 
maneuver as Fig. 9). In general, the histories march well. 
The computed time histories of q show a small bias from 
the measured data, but this did not affect the coefficient 
estimates or uncertainty. The computed time histories of 
N and V fit well with the measured data. The LSE and 
MLE algorithms were shown to satisfactorily match out- 
put time histories for well-conditioned maneuvers. 

Figure 12 presents the results of the subsonic maneuver 
that exhibited the smallest change in excess thrust (F ex ) 
from the nominal schedule at a flight condition of 0.85 
Mach and 35,000 ft. Time histories are presented for the 
following state parameters (F tx , M, ALT, a, 0 , q, N x , 

N ), and control parameters, (5 CJ , 8 4 ). Stabilized trim 
Hata were gathered from 0 to 8 sec. The steady-state 
value of N x was 0.031 g indicating that a bias was present 
or that the aircraft was slightly accelerating before the per- 
turbations (which did not affect the results). The nominal 
cowl deflection and stabilator trim position were 4.8° and 
0.62°, respectively. From 8 to 45 sec, the first trim of -3 
was sent to the cowl, causing the autopilot to drive die sta- 
bilator 0.23° trailing edge down to a new trim position of 
0.85°. This reduced N x and F tx by 0.006 g and 225 lbf, 
respectively, increasing the trim drag on the aircraft 
From 49 to 67 sec, the trim was removed and additional 


stabilized data were gathered. From 71 to 97 sec, the last 
trim of +3° was sent to the cowl, again causing the autopi- 
lot to drive the stabilator trailing edge up 0.37° to a new 
trim position of 0.25°. This trim increased N x and F tx by 
0.004 g and 90 lbf, respectively, decreasing the trim drag 
on the aircraft. 

This small change in N x was readily discernible in the 
time history and calculating a signal-to-noise ratio pro- 
vided an estimate of the smallest measured change m N x . 

This minimum was calculated to be 0.002 g, which 
equates to approximately 50 lbf of drag and coincides with 
resolution of the accelerometer. Note that the sinusoidal 
disturbance at a frequency of 0.32 Hz superimposed on the 
data was caused by an autopilot problem. Peak-to-peak 
changes of normal acceleration exceeded 0.40 g because 
of cyclical stabilator commands. As a result, it was drought 
at first that the data would be unusable. The postflight 
analysis, however, showed that both algorithms were 
robust to this type of disturbance and accurately estimated 
the gradients. For the LSE algorithm to be this robust was 
unexpected. 

Figure 13 is a summary plot of the three maneuvers per- 
formed during the investigation. The dimensional deriva- 
tive of longitudinal force coefficient with respect to cowl 
deflection is plotted as a function of Mach number. As 
was expected, the low dynamic pressure at the subsonic 
condition yields values close to zero, increasing as Mach 
number and dynamic pressure increase. Another area of 
interest that was not investigated is the mid-transomc 
regime, where difficulties in predicting an accurate cowl 
schedule could yield large gains over the baseline cowl 
schedule. There is overall agreement between the two 
analysis methods, MLE and LSE. The MLE approach has 
the added benefit of calculating Cram6r-Rao bounds or 
uncertainty estimates that would benefit a real-time control 

system. 

Figure 14 summarizes the demonstrated excess thrust 
benefits caused by incremental cowl rotation about its trim 
point Although the data for the three maneuvers flown do 
not represent the optimum cowl-stabilator configuration, 
they illu strate the performance gains from the first itera- 
tion to the optimum. The maneuver at Mach 1.25 and 

25.000 ft clearly shows that more than 800 lb of excess 
thrust were gained during the maneuver. The maneuver at 
Mach 1.6 and 35.000 ft was limited to gains of approxi- 
mately 250 lbf of excess thrust caused by shock ingestion 
constraints, but the gradient was large and similar to the 
maneuver at Mach 1.25. The maneuver at Mach 0.85 and 

35.000 ft was significant because it alone demonstrated a 
chiinge in sign of the gradient. Maximum excess thrust of 
100 lbf was estimated to occur at 8.5°, or T up from the 
9.5° pertuibation. 
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The success encountered in this feasibility study shows 
that either method could be used for real-time adaptive 
control of the cowl and encourages further research to 
app!y this approach to the engine itself. Given the similar 
results, faster computation, and easier implementation, the 
LSE method is more suitable for a real-time adaptive con- 
trol environment. Further studies also can determine 
whether the performance sensitivities produced by pertur- 
bations of the engine controls, compressor-inlet variable 
vanes, rear-compressor variable vanes, and nozzle area 
can be identified. 

Concluding Remarks 

The cowl was perturbed on the NASA F-15 highly inte- 
grated digital electronic control aircraft at three flight con- 
ditions to determine sensitivities of excess thrust to cowl 
deflections. The flight data were analyzed postflight using 
two estimation algorithms to compare accuracy, speed, 
and ease of implementation into a real-time flight control 
system. The maneuvers were flown at three conditions, 
one subsonic and two supersonic. 

The flight test results show that performance gradients 
always were identifiable with either method and that sub- 
stantial benefits are possible over the baseline cowl sched- 
ules, particularly at high Mach numbers. Both approaches 

the least squares estimator (LSE) and maximum likeli- 
hood estimator (MLE), produced similar results. 

The MLE approach yields more information, including 
Cramdr-Rao bounds, and contains a more complex three- 
degrees-of-freedom model that requires a more computa- 
tionally intensive iterative solution. The LSE approach 
contains a one-degree-of-freedom model that can intro- 
duce more inaccuracies into the estimate but can be solved 
directly. More expertise also is needed to initialize the 
MLE algorithm to ensure proper convergence, as opposed 
to the LSE method. 

Given the similar results, faster computation, and easier 
implementation, the LSE method appears more suitable 
for a real-time adaptive control environment. For adaptive 
implementation, only small perturbations that would not 
be noticed by either the pilot or passengers are required. 
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Appendix 

Ttie model sttucture used for the less! squarMflnilysK 
is described as the autoregressive exogenous (ARX) 
model where AR refers to the autoregressive part 
A(a)y U) and X refers to the extra input B ( q ) « (0 
called the exogenous variable. 1 Figure 1 5 is the oc 
diagram for the model. The general input-output relauon- 
shipis described by the foUowing linear difference equa- 
tion where n k is equal to the number of delays from input 
to output: 


y(t) + a x y(t- 1 ) + ••• + a ny( f 

b x u(t-n k ) +b 2 u(t-n k - D+ — 

+ b H u(t-n k - n b + 1 ) +<(0 ^ 9 ) 


where <f 1 is the delay operator, we see that the general 
predictor model 

y(!) = G{q,Q)n(t-n k ) +H(q,Q)e(t) (&) 


can be formed when 


G(o.e) = 


B(q) 


r 1 


H(q,Q) = 


(24) 


By definition, the noise characteristics of the - LSE are 
Gaussian Although it may not be intuitive that white 
noise is assumed to go through the denominator dynamics 
of the system before being added to the output, this model 
SL 7 ^importa* property to. makes i. a prune choree 
for many appUcations— the predictor defines a linear 
regression. This is shown below. 

The one-step-ahead prediction for equation (23) is given 


9 (,| 0 ) = /r‘ ( q , e) G ( q , 0)u(t- n k ) ^ 

+ [i-/r 1 («,e)]y (0 

Substituting equation (24) into (25) yields 

9(0 = B(q)u(t-n k ) + [1-A(<7)]y(0 


Now introduce the vector 


The estimated coefficients are contained in the following 
vector: 


0 = [a 






( 20 ) 


<p (0 = [-y(t-l).--y( f - rt a) 

u{t-n k )...u(t-n k -n b + 1 ) 1 


(27) 


Then equation (26) can be rewritten as 


If we introduce 

A(q) = \+a x q~ l + - + a H q 


901 0) = 0 T <p(O = <P r (06 ^ 28 - 

This predictor model is a linear regression composed of a 
(21) scalar product between the known data vector <P (0 and 

the parameter vector 0 . 


and 

B(q) = b 1 + M -1 + " , + i V 7 


( 22 ) 
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Fig. 3 The advanced airframe-propulsion technology system. 











Turbulence Noise 



Fig. 7 The MLE model. 
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Fig. 8 Test conditions for the performance sensitivity feasibility study. 
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Fig. 10 Measured and computed time histories of N using the LSE algorithm; Mach 1.25 at 25,000 ft. 
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Fig. 14 Summary of measurement-based performance optimization feasibility study cowl perturbation. 
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Fig. 15 The LSE autoregressive exogenous model. 
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